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ABSTRACT 

SGR 0501+4516 was discovered with the Swift satellite on 2008 August 22, after it emitted a 
series of very energetic bursts. Since then, the source was extensively monitored with Swift, the 
Rossi X-ray Timing Explorer (RXTE) and observed with Chandra and XMM-Newton, providing a 
wealth of information about its outburst behavior and burst induced changes of its persistent X-ray 
emission. Here we report the most accurate location of SGR 0501+4516 (with an accuracy of 0.11") 
derived with Chandra. Using the combined RXTE), Swift/X-r&y Telescope, Chandra and XMM- 
Newton observations we construct a phase connected timing solution with the longest time baseline 
(^240 days) to date for the source. We find that the pulse profile of the source is energy dependent 
and exhibits remarkable variations associated with the SGR 0501+4516 bursting activity. We also 
find significant spectral evolution (hardening) of the source persistent emission associated with bursts. 
Finally, we discuss the consequences of the SGR 0501+4516 proximity to the supernova remnant, SNR 
G160.9+2.6 (HB9). 

Subject headings: pulsars: individual (SGR 0501+4516 ) — X-rays: bursts 



1. INTRODUCTION 

Soft Gamma Repeaters (SGRs) are intriguing mani- 
festations of neutron stars. They are most often discov- 
ered in hard X-rays/soft 7-rays, when they (repeatedly) 
emit clusters of energetic bursts with peculiar proper- 
ties; these episodes are characterized by a plethora of 
'short bursts' that usually last only a fraction of a second 
but emit extremely large energies of 10 37 — 10 40 erg. On 
rare occasions, SGRs emit relatively more energetic, so 
called 'intermediate events', which are 1-2 orders of mag- 
nitude larger in energy, last longer (~l-4 s) and often ex- 
hibit an extended tail. On extremely rare instances (only 
three have been detected thus far), SGRs produce 'giant 
flares', which are even more energetic (10 44 — 10 47 erg) 
and have distinctive and similar morphologies: a very 
bright and spectrally hard initial short spike (< 0.5-s) 
that is followed by a longer (300—600 s) tail clearly mod- 
ulated at the spin frequenc y of the originating source (see 
iWoods k Thompsonll2006l , for a review). 

Along with their bursting behavior, SGRs and their 
close relatives, Anomalous X-ray Pulsars (AXPs) dis- 
play interesting persistent X-ray emission properties: 
they pulse with spin periods in the narrow range of 
2—12 s. They all slow down at rather large spin- 
down rates (~T0 -11 s/s), and their inferred dipole mag- 
netic fields are extreme ly strong (B-J3 ~ 10 14 — 10 15 G; 
iKouveliotou et al"1ll998T) . Moreover, their X-ray intensi- 
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ties vary, usually in the form of a rapid increase associ- 
ated with the onset of a major bursting episode, followed 
by a gradual decre a se lasting months to years (see e.g., 
I Woods et all 12001 lEsposito et all 1200a IWoods et all 
2004). 

SGRs and AXPs are now classified as magnetars - neu- 
tron stars that are powere d by their extremely stron g 
magnetic fields (B >10 14 G: lDuncan k Thompson! 1992f) . 
The magnetar model attributes the origin of the energetic 
SGR events to cracking of the solid neutron star crust, 
when strained by the drifting stro ng magnetic field: local 
fractu res leading to short bursts ((Thompson k Duncan! 
1995) and cracking at global scales giv ing rise to giant 
flares ([Thompson k Duncanl[l995ll200lH . The persistent 
emission of magnetars likely originates from the stellar 
surface heated by the decay of the strong magnetic field 
([Thompson k Duncadll996t lOzelll200ll : IHo k Laill20Qll ) 
and from resonant Compton scattering of soft photons by 
magnetospheri c currents driven by twists in the evolvin g 
magnetic field (|Thompson. Lvutikov. k Kulkarni r2002). 

SGR 0501+4516 was discovered on 2008 August 22, 
when four short and soft bursts from the source trig- 
gered the Burst Alert Telescope (BAT) onboard NASA's 
Swift satellite (B arthelmv et al.l 120081 ). To confirm its 
magnetar nature, we initiated our Target of Opportu- 
nity (ToO) observations with the Rossi X-ray Timing Ex- 
plorer (RXTE). The first iJYTS/Proportional Counter 
Array (PCA) observation took place on August 22 for 
a net exposure of 484 s. Timing analysis of this short 
stretch of dat a revealed coherent pulsations with a pe- 
riod of 5.76 s (Gogiis et al. 2008). Subsequent observa- 
tions of the new source with RXTE an d the Swift/X-ray 
Telescope (XRT. iBurrows et al.|[2005h allowed us to de- 
termine the spindown rate of SGR 0501+45 16, and there- 
fore, to firmly establi sh it as a magnetar (|Woods et al.l 
2008; IRea et al.ll2009T) . 

Following the onset of its outburst episode, the source 
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flux decaye d exponentially with an e-folding time of 
23.8 days (|Rea et al.l I2009D . Using a subset of the 
Swift/XRT observations we also employ here and in 
addition their X MM-Newton and Suzaku observations, 
IRea et ahl ()2009D determined a pulse ephemeris, which 
showed evidence of strong negative P. The latter 
was interpreted as the recovery to a secular spindown, 
which might have incre a sed in conjunction with the out- 
burst onset. IRea et aTl (|2009( ) also analysed an archival 
ROSAT observation of the source and found that the 
X-ray flux of SGR 0501+4516 in 1992 was about 80 
times lower. Recently, using Suzaku observations of 
SGR 0501+4516 obtained four days after the outburst 
onset, lEnoto et al.l (|2010f ) reported that the persistent 
emission of the source extends up to about 70 keV. 
Finally, besides Swift/BAT, SGR 0501+4516 also trig- 
gered the Gamma ray Burst Monitor (GBM) onboard 
the Fermi Gamma-ray Space Telescope (FGST) 26 times. 
Both instruments recorded burst activity from the source 
until about 2008 September 3. Detailed analysis of the 
GBM data is currently underway (L. Lin et al., in prepa- 
ration) . 

We present here spatial, and long term temporal and 
spectral characteristics of the persistent X-ray emission 
of SGR 0501+4516 using wideband X-ray observations 
during three states of the source: (i) its burst active 
episode, (ii) its transition to quiescence and finally, (iii) 
its quiescent phase. In Section 2, we describe the ob- 
servations and data used in our study. In Section 3, we 
present the results of our detailed analyses. We discuss 
the interpretation of our results in Section 4. 

2. OBSERVATIONS 
2.1. RXTE 

We monitored the source with RXTE until 2008 Au- 
gust 30 in densely spaced pointings (consecutive inter- 
vals ranging between 0.1 — 0.9 days) followed by longer 
intervals as the source became fainter (spaced 2—7 days 
apart). Exposure times of these observations varied be- 
tween ~1 ks (in two cases) and ~11.6 ks (in one case), 
while the great majority lasted ~2.5 ks. Overall, we ac- 
quired a total exposure of 82.5 ks in 29 RXTE pointings 
from 2008 August 22 (MJD 54700) to October 14 (MJD 
54753). We used data in collected with the RXTE /PC A 
only. The number of operating proportional counter 
units (PCUs) varied between 1 and 3, with a median 
value of 2 during our pointings. For each observation, 
we ap plied our burst search algorithm (see iGogiis et al.l 
(2000) for its methodology) to the lightcurve with 0.125 
s time bins to identify short events. We then filtered out 
the times of short bursts identified to obtain a burst-free 
event list for timing studies. We finally converted all 
event arrival times to the Solar system barycenter. 

2.2. Swift 

SGR 0501+4516 was monitored with the Sunft/XRT 
starting 2008 August 26 (MJD 54704) until 2009 April 
19 (MJD 54940) for a total exposure of 436 ks, provid- 
ing the longest baseline coverage of the source to date. In 
Figure[T]we present the long term X-ray light curve varia- 
tions of SGR 0501+4516 as seen in the 0.3-10 keV range 
with Swift/ XRT. The solid line is obtained by fitting an 
exponential function, whose time reference is the time 
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Fig. 1— Swift/XRT X-ray (0.3-10 keV) light curve of 
SGR 0501+4516. The solid line is the best fit exponential function. 
The dates shown on top are in YY/MM/DD format. 

of the first burst observed on August 22 (i.e., 54700.529 
MJD), that yields an e-folding time of 27.9+2.5 dayfl 
This value i s notably consist ent (within 1.6cr) with that 
reported by IRea et al.1 (|2009f) using the first 160 days of 
data coverage of the source. 

We used all XRT observations that were performed in 
Windowed Timing (WT) mode. For each pointing we 
extracted source events from a 15"segment of the one 
dimensional WT image, centered at the pixel with the 
highest rate. The background events were selected from 
source free intervals on both sides of the source extraction 
interval. We then generated and filtered the light curves 
of each XRT pointing to ensure that there are no short 
SGR bursts in the persistent emission data included in 
our timing analysis. Finally, we transformed the arrival 
time of each event to the Solar system barycenter for 
timing studies. 

We also employed additional Swift/XRT observations 
in Photon Counting (PC) mode that took place around 
2009 October 9, December 6 and 2010 February 20 for 
13.2 ks, 13.4 ks and 11.6 ks, respectively. We performed 
spectral analysis with these data sets and derived the 
source flux at late stages of the outburst decay. 

2.3. Chandra X-ray Observatory 

We initiated our target of opportunity observations 
with the CTitmrfra/Advanced C CD Imaging Spectrome- 
ter (ACIS: iGarmire et al.ll2003T l in Continuous Clocking 
(CC) mode on 2008 August 26 (MJD 54704) for an ef- 
fective exposure of 36.5 ks. The source was relatively 
bright at the time of the ACIS observation resulting 
in a (background subtracted) count rate of 4.42+0.01 
counts/s (0.3-8.5 keV), which is still below the level that 

9 Note that uncertainties reported throughout the paper are lcr 
unless otherwise indicated. 
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would cause pile-up. We extracted source events from an 
8"segment of the 1-dimensional CC mode image centered 
at the pixel with the highest counts; we extracted back- 
ground events from a 20"long segment with a 10"buffcr 
on both sides of the source extraction. We then applied 
barycentric correction to all event arrival times. 

We also observed SGR 0501+45 16 with the Chan - 
dra/ High Resolution Camera f#RC; lMurrav et al.ll2000h 
in Imaging mode for 10 ks on 2008 September 25 (MJD 
54734). We used the superb angular resolution of the 
HRC to derive the precise X-ray position of the source. 
The X-ray intensity of the source was considerably lower 
at the time of the HRC observations with a (background 
subtracted) count rate of 0.505±0.007 counts/s in the 
0.3—8.5 keV range. 

2.4. XMM-Newton 

We observed SGR 0501+4516 with XMM-Newton on 
2008 August 23 (MJD 54701) for 48.5 ks. Here, we used 
data collec ted with the European P hoton Imaging Cam- 
era (EPIC; Pfeffermann "eTailll999h PN camera operated 
in Imaging with Prime Small Window mode. We se- 
lected source events from a circular region with a radius 
of 20"centered at the source and background events from 
a source- free region on the same chip. The background 
subtracted source count rate is 5.62+0.01 counts/s in the 
0.3—8.5 keV range. Inspection of the source light curve 
revealed many short SGR events, as expected since the 
observation took place during the burst active episode 
of the source. We adopted our burst finding algorithm 
to search in the 0.3-10 keV band of the EPIC-PN data 
and identified 642 short bursts over a wide range of X- 
ray intensities. After removing the time intervals corre- 
sponding to all bursts, the exposure time was reduced to 
47.5 ks. Finally we applied a barycentric correction to 
the event arrival times for our timing studies. 

3. DATA ANALYSES AND RESULTS 
3.1. Source Location 

We generated an image of the entire Chandra/ HRC- 
I field in the 0.5—7.0 keV range and used wavde- 
teclH to identify all X-ray sources in our frame. We 
found two point sources: SGR 0501+4516 at the aim 
point of the image and a previously uncatalogued X- 
ray source (CXOU J050111. 3 + 451525) that is separated 
by 1.4' from the SGR location. Their locations are 
RA = 05 ft 01 m 06 s .76, Dec = +45°16'33.92" (J2000), and 
RA = 05' l 01 m ll s .33, Dec = +45°15'25.28" (J2000) for 
SGR 0501+4516 and CXOU J050111.3 + 451525, respec- 
tively. 

The statistical uncertainty of the X-ray location of the 
latter source is 0.09" (22 source counts). To derive a 
precise positional uncertainty, we searched the 2MASS 
catalogue and found an IR source (2MASS 05011132 + 
4515252) at RA = 05 /l 01 m ll s .33, Dec = +45°15'25.25" 
with an uncertainty of 0.06". This location is consistent 
within the uncertainties, with the Chandra location of 
CXOU J050111.3 + 451525, so we concluded that the 
latter is very likely the X-ray counterpart to 2MASS 
05011132 + 4515252. Given the excellent alignment be- 
tween these two localizations, we estimate the absolute 

11 a CIAO tool, http://cxc.harvard.edu/ciao/ 



positional uncertainty of SGR 0501+4516 to be 0.11" 
(la). 

3.2. Pulse Properties 
3.2.1. Pulse Timing Analysis 

For this analysis we used the puls e timin g technique 
described in detail in iWoods et al.l (|2002j ). In sum- 
mary, we use an epoch folding algorithm to determine 
the pulse frequency and its higher order time deriva- 
tives as follows. We first group all available observa- 
tions (RXTE/PCA, Swift/ XRT, Chandra/ ACIS, XMM- 
Newton/EPIC-PN) into discrete time segments. We then 
generate a (high signal to noise ratio) pulse profile tem- 
plate using the data in the 2 — 10 keV energy band by 
accumulating multiple, closely spaced pointings. Next we 
cross-correlate the pulse profile obtained by each group 
of pointings with the template profile and measure the 
phase drifts with respect to the template (see Figure [2j 
top panel) . We obtained a good fit to the phase residuals 
by modeling with a 5th order polynomial (x 2 /degrees of 
freedom = 139.5/110). In Table [1] we provide the best 
fit parameters of the spin ephemeris of SGR 0501+4516. 
The fit residuals are presented in the middle panel of 
Figure [2j We also employed a quadratic spline fit to all 
available phase residuals in seven segments. In Table [2] 
we list the spin frequencies and their time derivatives as 
determined with the quadratic spline method; the bot- 
tom panel of Figure [5] shows the phase residuals obtained 
with the spline fit. It is important to note that in the 
pulse timing technique we employ, a systematic pulse 
profile change in which the profile drifts in phase cannot 
be distinguished from phase drift due to timing noise. 

TABLE 1 

Pulse ephemeris of SGR 0501+4516 obtained using the 
combined RXTE/PCA, Swift/XRT, CXO/ACIS-S and 
XMM-Newton/ 'EPIC-PN observations 



Parameter 


Value 


Range (MJD) 


54700.794 - 54940.953 


Epoch (MJD) 


54750.0 


v (Hz) 


0.173547943(1) 


v (10~ 13 Hz s- 1 ) 


-1.752(8) 


v (10~ 21 Hz s~ 2 ) 


6.9(3) 


„(3) (io-27 Hz s" 3 ) 


-2.7(3) 


„(4) ( 10 -34 Hz s -4) 


3.5(6) 



TABLE 2 

Pulse frequency and frequency derivative of 
SGR 0501+4516 obtained with a quadratic spline fit to 
the RXTE/PCA, Swift/XRT, Chandra/ 'ACIS-S and 
XMM-Newton/PN DATA 



Epoch 


Range 


V 


v 


(MJD) 


(MJD) 


(Hz) 


(10- 13 Hz s" 1 ) 


54707.5 


54700.0 - 54715.0 


0.173548659(2) 


-2.03(8) 


54722.5 


54715.0 - 54730.0 


0.173548384(3) 


-2.21(8) 


54740.0 


54730.0 - 54750.0 


0.173548100(3) 


-1.62(8) 


54765.0 


54750.0 - 54780.0 


0.173547721(4) 


-1.85(6) 


54805.0 


54780.0 - 54830.0 


0.173547146(3) 


-1.55(4) 


54855.0 


54830.0 - 54880.0 


0.173546394(4) 


-1.93(4) 


54915.0 


54880.0 - 54950.0 


0.17354545(2) 


-1.75(6) 
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Fig. 2. — Top panel Plot of pulse phase drifts of each 
SGR 0501+4516 observation with respect to the pulse template 
(see also the text); Middle Panel the fit residuals when the pulse 
phase evolution is modeled with a fifth order polynomial; Bottom 
Panel the fit residuals when a quadratic spline fit is applied. 

3.2.2. Pulse Profile Evolution 

To investigate the pulse morphology and its evolution 
in energy and time, we generated pulse profiles in the 
0.3-2 keV, 2-4.5 keV, and 4.5-8.5 keV energy inter- 
vals for the Chandra/ACIS-S, XMM-Newton/EPIC-PN 
and Swift/XRT data, and in the 2-4.5 keV, 4.5-8.5 keV, 
8.5-14 keV and 14-40 keV intervals for the RXTE/PCA 
data. We used the phase connected spin ephemeris pre- 
sented in Table [T] for determining the phases of all event 
arrival times in each specified energy interval. 

Figure [3] shows the pulse profile evolution with en- 
ergy as observed with XMM-Newton/EPIC-PN on the 
day after the onset of activation and during an episode 
of continued bursting. We find that the pulse profile 
in the 0.3—2 keV band is dominated by a broad peak 
(0-0.25-0.75) and a broad valley (0-0.75-1.25), with 
overlying sub-pulses. The most significant sub-pulse ap- 
pears at 0—0.34. The significance of these sub-pulses 
decreases at higher energies and the profile becomes al- 
most sinusoidal in the 4.5—8.5 keV range. The bottom 
panel in Figure [3] demonstrates that the lowest energy 
interval emission is the most dominant component of the 
overall profile. 

The RMS pulsed fractiorQ is energy dependent 
and increases with energy: (24.2±0.2)%, (30.6±0.3)%, 
(32.9±0.8)% in 0.3-2 keV, 2-4.5 keV, 4.5-8.5 keV, re- 
spectively. The RMS pulsed fraction in the entire energy 
range (0.3-8.5 keV) is (27.1±0.2)%. 

In Figure 01 we show the pulse profiles of 
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Fig. 3.— XMM-Newton/EPIC-PN pulse profiles of 
SGR 0501+4516 with their RMS pulsed fraction values in 
various energy ranges as indicated on each panel. 

SGR 0501+4516 seen with Chandra/ACIS-S near the 
end of the burst active episode of the source. We find 
that the morphology of these profiles closely resemble in 
all energy ranges those obtained with the XMM-Newton, 
despite important differences in the source state: (i) the 
intensities (count rates) in all energy ranges are lower, (ii) 
the sub-pulse that peaks around 0—0.03 becomes more 
prominent, and (iii) some pulsed fractions are signifi- 
cantly different. Most prominently, the Chandra 0.3—2 
keV RMS pulse fraction is (21.6±0.3)%, which is 7.2ct 
lower than that seen with XMM-Newton. The 2-4.5 keV 
range pulse fraction is (29.7±0.4)% (consistent with the 
XMM-Newton data within error). Finally the 4.5-8.5 
keV band has an RMS pulsed fraction of (38±1)%, which 
is 4cr higher than that in the XMM-Newton data. The 
energy-averaged RMS pulsed fraction (0.3—8.5 keV) is 
(25.2±0.3)%. 

We present in Figure M the RXTE/PCA energy- 
resolved pulse profiles and their evolution in time. For 
this analysis, we have split the PGA observations into 
4 groups of nearly equal exposures and obtained pulse 
profiles for each in the instrument energy range of 2 — 40 
keV. As the PCA is not an imaging instrument, its back- 
ground determination relies on model based estimates 
and most likely does not reflect the true background; for 
that reason, we plot the pulse profiles in arbitrary units. 

The leftmost column of Figure [5] corresponds to pulse 
profiles obtained during the most burst active episode 
of the source (—75% of bursts were emitted in the two 
days following the episode onset). The pulse profile in 
the 2—4.5 keV range is characterized by a single, smooth 
pulse with a large duty cycle (—60%). In the 4.5—8.5 
keV band, the pulse structure is still broad and peaks 
around 0—0.5, as in the lower energy interval, but it is 
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Fig. 4.— Chandra/ ACIS-S pulse profiles of SGR 0501+4516 with 
their RMS pulsed fraction values in various energy ranges as indi- 
cated on each panel. 

asymmetric. Note that there is also evidence of weaker 
structures around <£~1.0. The pulse profile in the 8.5 — 14 
keV band is multi-peaked: the structure around </>~1.0 
becomes more prominent and almost as significant as the 
other pulse in the profile at about half the spin cycle. 
However, only the structure around 0~1.O remains in 
the 14—40 keV energy band. 

In the time interval from 2009 August 25 to 30, (the 
second column from the left in Figure [5]) , the source was 
at a very low level of bursting activity. Here we see again 
the structure at around 0^1.0 and the sub-pulse near 
</>~0.34 that was also seen in the 2—4.5 keV band pro- 
files of XMM-Newton and Chandra. In the higher energy 
bands there are no remarkable changes in pulse shapes 
compared to the earlier interval, except for the fact that 
the structure near </>~1.0 is now less significant. 

The last short burst from SGR 0501+4516 was ob- 
served on 2009 September 3, therefore, the profiles cor- 
responding to the time intervals from September 6 to 
24, (third column from the left in Figure [5]), and from 
September 27 to October 14, (the rightmost column in 
Figure EJ, are characteristic of the source's descent to 
the quiescent phase. In the 2—4.5 keV band of the third 
column we find a sub-pulse around <f>~0.7 that was sig- 
nificantly seen in the 0.3—2 keV Chandra profile. The 
structures peaking at 4>~0.5 and </>~1.0 are clearly seen 
in the 4.5—8.5 keV band, while they are marginally seen 
in the 8.5—14 keV range. The pulse profile of the high- 
est energy interval only shows marginal evidence of the 
structure near </>~1.0. The 2—4.5 keV pulse profile of the 
rightmost column contains a single structure with a large 
duty cycle, broadly resembling profiles in the same band 
of earlier intervals. We find that the duty cycle of this 
pulse is significantly lower (^0.35 of the spin phase) in 



the 4.5—8.5 keV energy band. The pulse profiles of the 
higher energy bands of this episode are consistent with 
random fluctuations. 

Finally, we also studied the SGR 0501+4516 en- 
ergy and time resolved pulse profile evolution using the 
Swift/XRT observations which span a much longer time. 
We grouped the XRT pointings into four segments, each 
of which is comprised of nearly equal exposure time. The 
first interval corresponds to the burst active episode of 
the source, (the leftmost column in Figure [6]), the second 
interval starts immediately after bursting ceased, (the 
second column from left in Figure |6]), while the third and 
fourth intervals correspond to the quiescent state (long 
after the burst active phase; the observation end date is 
2009 February 4). Unlike the RXTE/PCA pulse profiles 
in Figure we plot the Swift/XRT profiles in real units 
since proper background subtraction is possible with the 
XRT data. We find that the sub-pulses seen in the low 
energy bands of the XMM-Newton and Chandra obser- 
vations persist during the burst active episode and in the 
interval immediately after, but disappear during quies- 
cence. Interestingly, the sub-pulse near </>~0.7 becomes 
remarkably prominent in the 0.3—2 keV and 2—4.5 keV 
ranges within days as the bursting stops. The pulses pro- 
files of the interval of 40—65 days after the outburst on- 
set reveal that the pulsed amplitudes in all energy ranges 
were significantly reduced (the third column from left in 
Figure |6|) and sub-pulses seen in the lower energies have 
disappeared. During the quiescent phase (the rightmost 
column in Figure |6|) , we find further gradual decline in 
pulsed intensity, while similar to the previous time in- 
terval, we see only one significant pulse with duty cycle 
decreasing with energy. 

We plot in Figure Q the energy resolved variations of 
the RMS pulsed fractions during the burst active episode 
and immediately afterwards in the 0.3—2 keV, 2—4.5 keV 
and 4.5—8.5 keV intervals with XMM-Newton, Chandra 
and Swift/XRT. In the 0.3-4.5 keV band the pulsed frac- 
tion drops to its minimum value during ongoing source 
bursting. The lowest energy band RMS increases with 
time, while in the other two bands the deviation from a 
constant trend is less than 2.4<r. 

3.2.3. Pulsed Intensity Properties 

To investigate the pulsed intensity variations of 
SGR 0501+4516 in time and energy, we employed the 
RMS pulsed count rat^3, which is not affected by back- 
ground, therefore provides a measure of the pulsed inten- 
sity. We used the PC A observations since its exclusive 
coverage starts from soon after the onset of the source 
activation, continues throughout the active episode and 
extends well into the burst quiescent phase, while the 
Swift/XRT observations start after the end of the major 
bursting activity. We present in Figure [8] the evolution 
of the RMS pulsed intensity in three energy intervals. 
We find that the pulsed intensity in the 2—4.5 keV band 
gradually declines while the source is actively emitting 
bursts, while the pulsed intensities in the 4.5—8.5 keV 
and 8.5—14 keV bands increase at the same time. The 

15 PCR rms = i (Eili(Ri -Rave) 2 - ARfY with propa- 
gated errors 5PCR rms = N PC 1 Rrm3 (e£=i [(Ri " Ravo)AR;] 2 ) 3 
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Fig. 5. — Pulse profile history of SGR 0501+4516 in the energy intervals shown on the right of each horizontal panel, obtained with 
RXTE/PCA. The labels on top of columns are the time ranges (in MM/DD format) within which the data used to construct these profiles 
were accumulated. 
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Fig. 6. — Pulse profile history of SGR 0501+4516 in the energy intervals shown on the right of each horizontal panel, obtained with 
Swift/XRT. The labels on top of columns are the time ranges (in MM/DD format) within which the data used to construct these profiles 
were accumulated. Note that the year changes to 2009 at the rightmost column. 



pulsed intensities remained constant in all energy bands 
for about a week following the active bursting phase, and 
then declined. The steepness of the decline is inversely 
proportional to the energy range. 

In Figure [5] we show the evolution of the pulsed hard- 
ness ratios, which are the ratios of the pulsed intensity 
in the two higher energy bands to that in the 2—4.5 keV 
band. We find that both ratios increase rapidly in the 
first 1.5 days into the outburst, continue to increase at 
a slower pace until the end of the burst active phase, re- 
main constant for about a week, and decline afterwards. 

3.3. Spectral Properties 

We analyzed the X-ray spectra of SGR 0501+4516 us- 
ing the XMM-Newton, Chandra/ AC IS- S and the 
Swift/XRT observations during and immediately after 



the burst active episode of the source. With respect 
to the longest extending observations of the Swift/XRT, 
we generated 7 combined spectra using multiple point- 
ings to provide statistically significant spectral parame- 
ters and close in time, to avoid blending different spec- 
tral shapes (in case of spectral evolution). We fit all 
spectra simultaneously with an absorbed blackbody plus 
power-law model by linking the interstellar absorption 
parameter in the 0.5—10 keV energy range. We obtained 
an acceptable fit to all spectra (x 2 /degrees of freedom 
= 5086.2/4470 = 1.14) with a hydrogen column density, 
iV H =(0.95±0.08)xl0 22 cm" 2 . In Figure [TO] we present 
the evolution in time of our best fit blackbody temper- 
ature and power-law index. We find that the blackbody 
temperature declined very rapidly while the source was 
burst active, increased to a maximum value of kT = 0.77 
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Fig. 7. — Energy resolved RMS pulsed fraction variations 
of SGR 0501+4516 as obtained with the data of XMM- 
Newton/EPIC-PN, Chandra/ ACIS-S and Swift/XRT. The vertical 
dashed lines indicate the interval during which multiple short SGR 
bursts were detected. Letters at the bottom of the Figure indicate 
the instrument with which the plotted data above were obtained; 
X: XMM-Newton, C: Chandra and S: Swift/XRT. 



Time (days since the outburst onset) 

Fig. 8. — Time history of pulsed count rates in the 2—4.5 keV 
(filled circles), 4.5—8.5 keV (squares) and 8.5—14 keV (diamonds). 
Note that the tick mark values on the left axis correspond to the 
2—4.5 keV band data, while the values on the right axis to the 
higher two energy bands. The vertical dashed line indicates the 
end of the active bursting episode of the source. 



keV soon after the intense bursting ceased, and then 
gradually decreased. The power law component, on the 
other hand, became less steep during the burst active 
phase, and then steeper as the source went into burst 
quiescence. 

Motivated by the rapid spectral variations (as well 
as by the change in RMS pulsed fraction) seen in the 
Chandra/ ACIS-S observations, we performed a phase- 
resolved spectral analysis of the data. As a reminder, 
the Chandra data cover the fourth day after the onset of 
the source burst activity. We divided the spin phase 
into 10 bins and extracted spectra for each bin. We 
then simultaneously fit all 10 spectra with the attenu- 
ated blackbody plus power-law model, again by linking 
the interstellar absorption. We obtained a good fit to all 
spectra (x 2 /degrees of freedom = 1909.7/1785 = 1.06) 
with N H =(0.98±0.03)xl0 22 cm" 2 , consistent with that 
obtained for the integrated spectra. In Figure [IIJ we 
present the best fit phase-resolved blackbody tempera- 
tures (left) and power law indices (right). We find that 
the blackbody temperature remains constant within er- 
rors over the whole spin phase, while the power law com- 
ponent is strongly intensity-dependent: the spectrum be- 
comes softer during the phase minimum and harder dur- 
ing the phase maximum. 

Finally, to evaluate the variations in spectral shape 
and source flux as the source entered the quiescent 
phase, we generated 3 Swift/XRT spectra obtained with 
about 12—13 ks of exposures that took place long af- 
ter the outburst episode (2009 October, December and 
2010 February). We fit all three spectra simulta- 



neously with the absorbed blackbody plus power-law 
model and with the hydrogen column density linked. 
We find that all three power-law indices are consis- 
tent with one another, therefore, we linked them as 
well to better constrain the blackbody component. We 
obtained a very good fit (x 2 /degrees of freedom = 
123.8/114 = 1.09) with a iV H =(1.19±0.22) x 10 22 cur 2 , 
a common power-low index of 5.04±0.21 and black- 
body temperatures of 0.78±0.07 keV, 0.59±0.07 keV and 
0.56±0.06 keV, for the above mentioned three epochs, 
respectively. The corresponding 1 — 10 keV unabsorbed 
flux values were (7.0±0.3) x 10~ 12 , (6.4±0.2)xl0~ 12 and 
(5.8±0.2)xl0~ 12 erg cm" 2 s _1 . 

4. DISCUSSION 
4.1. Persistent emission properties 

We determined the pulse ephemeris of 
SGR 0501+4516 using the longest time baseline 
available (^240 days). We find that the spin ephemeris 
of the source can be described well with a fifth order 
polynomial (sec Table 1) or by employing a quadratic 
spline method including only the spin frequency and 
frequency derivative (Table 2). We, therefore, argue 
that the recent rep ort of P being t he recovery of the 
secular spin-down (|Rea et al.l [2009) is not the only 
viable possibility. The two best studied SGR sources 
(SGR 1806-20 and S GR 1900 + 14) are known to 
exhibit torque noise dWoods et all (20021); similarly, 
AXP1E 1048 . 1-593 7 (jGavriil fc Kaspil l2004). Recently, 
iBeloborodovl (|2009f ) suggested that untwisting of the 
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Fig. 9. — Time history of hardness ratios of RMS pulsed inten- 
sities in the 4.5—8.5 keV to that in the 2—4.5 keV (filled circles) 
and in the 8.5—14 keV to that in the 2—4.5 keV (diamonds). The 
vertical dashed line indicates the end of active bursting episode of 
the source. 



magnetic field through the magnetosphere due to sudden 
crustal motion may account for the variable or noisy 
spin-down behavior of these sources. The variations 
seen in SGR 0501+4516 may then well be due to a noisy 
spin-down of the source. 

We estimate the inferred dipole magnetic field strength 
of SGR 0501+4516 as B d 2x10" G, assuming that 
the neutron star slows down via magnetic braking. This 
value is very similar to the inferred field strength of 
the r ecently discovered SGR J1833— 0832 (Gogi is et al.l 
2010) and lies near the lower end of the inferred mag- 
netic fields of SGRs. The characteristic age (P/2P) of 
SGR 0501+4516 is estimated to be 15.7 kyr. 

We find the energy and time dependence of the pulse 
profiles of SGR 0501+4516, during its burst active 
episode and soon afterwards, intriguing. During the 
burst active phase, the ~4.5 keV pulse profiles are 
smooth while those in the lower energy bands show sig- 
nificant subpulses. In particular, there are at least four 
subpulses between 0.3—2 keV, indicating that there are 
significant intensity variations over the spin phase. The 
most significant subpulse peaks near 0^0.34 and lasts 
about 0.72 s, which corresponds to ~l/8 of the spin pe- 
riod. 

The energy dependence of the pulse profile extends to 
the higher energies; the pulsations are clearly seen up to 
40 keV following the onset of the outburst. The signifi- 
cant peak in the 14—40 keV band near 0~1.O is weaker in 
the 8.5—14 keV range and not seen in the lower energy 
bands. This hard component fades a way within a few 
weeks after the burst activity onset. iRea et ail (|2009f ) 



Time (days since the first burst) 

Fig. 10. — Evolution of the X-ray spectral parameters of 
SGR 0501+4516: blackbody temperature (diamonds) and power- 
law index (squares). The vertical dashed lines indicate the interval 
during which intense short SGR bursts were detected. Letters at 
the bottom of the Figure indicate the instrument with which the 
plotted data above were obtained; X: XMM-Newton, C: Chandra 
and S: Swift/XRT. 
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Fig. 11. — Phase-resolved spectral parameters of 
SGR 0501+4516 measured using Chandra/ ACIS-S observa- 
tions. Left Panel: blackbody temperatures; Right Panel: power 
law indices. We overplot the pulse profile (histogram with units on 
the right) binned into phase steps of the phase resolved spectral 
extraction. 



and lEnoto et al.l (|2010l) deduces the same decay behav- 
ior using INTEGRAL and Suzaku/HXD observations, 
respectively. A similar trend is also seen in the post- 
outburst pulse profiles of SGR J0418+5716 (P. M. Woods 
et al. 2010, in preparation). This fact may (indirectly) 
indicate the existence of heating of the neutron star crust 
as a result of crustal fracturing (|Thompson fc Duncan! 
[19961) . 

We find that the pulsed fraction significantly increases 
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with energy after a large drop in the 2—4.5 keV band near 
the end of the bursting episode. This trend is noteworthy 
but should be taken with caution since inter-instrument 
calibration may play a role as we discuss below. 

SGR 0501+4516 went through significant spectral vari- 
ations both during the bursting phase and throughout 
the quiescent phase, when no more SGR events were seen. 
We find that the spectrum of th e source was har dening 
during the active period (see also lRea et al.l ((2009)). Ac- 
cording to the magnetar scenario, this is a natural out- 
come of surface fracturing at a wide range of scales: suffi- 
ciently large scale crus tal fracturing would lead to highly 
energetic short bursts (|Thompson fc Dunc an 1996) while 
smaller scale fracturing may con tribute to the persis- 
tent X-ray emission of the source (|Thompson fe; Duncan! 
1995). During the post-outburst episode, the spectrum 
of SGR 0501+4516 remained constant for about a week 
and softened thereafter. This is also expected in the 
magnetar model a s the heated stellar crust cools down 
(|Guver et al.ll2007ft . 

The significant phase-resolved spectral variations near 
the end of the burst active phase seen with Chandra are 
very interesting. We showed that the blackbody tem- 
perature remains constant over the spin phase, while the 
power law component gets shallower (harder) with in- 
creasing X-ray intensity. Since the peak interval of the 
pulse phase would dominate the phase-averaged spec- 
trum, it is possible to 'hide' a shallow power law trend. 
On the other hand, the sharp drop in the blackbody tem- 
perature during the burst active phase and then its rapid 
rise in a day is indeed puzzling. Since this behavior was 
only seen with the Chandra data, there is no way to con- 
firm these results with other contemporaneous observa- 
tions. 

Finally, we find that the flux in the 1 — 10 keV range 
measured with Swift/XRT between 412 and 546 days af- 
ter the onset of the 2008 outburst, follows a marginal de- 
cline from 7xl0~ 12 to about 6xl0~ 12 erg cm" 2 s _1 . It 
is interesting to note that the source flux in the same en- 
ergy range at about 100 days after the onset was around 
6— 7x l0~ 12 erg cm -2 s _1 (see Figure 11 of iRea et al.1 
12001 . Therefore, the source flux likely remained con- 
stant over a year following a steady decrease over the 
first 100 days after activation. 

4.2. Archival outburst episodes of SGR 0501+4516 

On 1993 July 25, the Burst And Transient Source Ex- 
periment (BATSE) onboard the Compton Gamma Ray 
Observatory (CGRO) triggered on two peculiar events: 
at 17:45:21 UT (BATSE trigger number 2463) and at 
23:32:37 UT (BATSE trigger number 2464). In Fig- 
ure [121 we present BATSE Large Area Detector (LAD) 
light curves of these two events in the four LAD dis- 
criminator energy intervals. These two bursts were both 
short (T 90 = 0.064 s and 0.048 s, respectively), had 
soft spectra and their locations were consistent with 
one another (RA, Dec= 05 h 09 m 12 s , 43°48' 00", er- 
ror radius of 5.7°, and 04 h 42 m 00 s , 44° 12' 00", error 
radius of 6.2°, respectively). Note that the Chandra 
position of SGR 0501+4516 falls well within the over- 
lap of these two error circles. Given the soft nature 
of these events, their short durations and the fact that 
they originated from the same region of the sky within 
about 5.8 hours apart, it is strongly suggestive that these 
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Fig. 12.— Light curves of BATSE Trigger number 2463 (left) 
and 2464 (right) in the four LAD discriminator energy intervals 
which are indicated on the right. 



events are from SGR 0501+4516 . Our search for un- 
triggered events using the c ontinuous BATSE DISCLA 
data (see lGogus et al.lll999l for the methodology) for 10 
days starting on 1993 July 20, resulted in no additional 
events from the source. Therefore, we conclude that 
SGR 0501+4516 was indeed burst active back in July 
1993 for a very brief period. 

Furthermore, the Gamma Ray Burst Monitor (GRBM) 
onboard the BeppoSAX satellite recorded a short burst 
(T 90 = 0.77 s) on 2000 October 11 at 10:06:49 UT (GRB 
001011A) from RA = 04 h 28 m , Dec = 53°(with a po- 
sitional uncertain ty of 34°) that inclu des the position of 
SGR 0501+4516 (jFrontera et al.|[2009h . Even though the 
spectrum of the GRBM burst is not very soft, there is an- 
other short duration event from the same direction about 
95 s after the GRBM trigger. It is, therefore, suggestive 
that the source of these bursts is SGR 0501+4516 since 
they were positionally coincident, they were short and 
they repeated. Note, however, that the large error cir- 
cle of GRBM also included another recently discovered 
magn etar source, SGR J0418+5729 (|van der Horst et all 
12010) . 

With the definite active episodes in 1993 and 2008, 
and a third possible one in 2000 we can infer that 
SGR 0501+4516 is going through burst active phases 
on a timescale of ^10 years. This timescale is 
similar to that of SGR 1627—41 which was discov- 
ered in 1998 during about a we e k-long burst active 
episode (|Kouveliotou et al.l 119981: iWoods et ail 119991) 
and exhibited a burst active epoch again in 2008 
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(|Esposito et al.ll2008f) . The two recently discovered SGRs 
(SGRJ0418+5729 and SGR J1833-0832) have emitted 
only a few bursts during their very brief activ e episodes 
([van der Horst et al.ll2010t iGogiis et all l2010f ) and their 
outburst recurrence behaviors may characteristically be 
different from SGR 0501+4516 and SGR 1627-41. ' 

4.3. On the association with Supernova Remnant, 
G 160. 9+2.6 

iGaensler fe Chatteriee I (|2008f ) comp ared the posit ion 
of SGR 0501+4516 with the locations (|Green 1120091) of 
known Galactic supernova remnants (SNRs), and re- 
ported that the magnetar sits just outside the south- 
eastern rim o f the well-known SNR G160.9+2.6 (also 
know n as HB9 Han burv fe Hazard~| [T953; Lea hy fe Tian I 
[2001 . 

All of the four initially identified SGRs (1806-20, 0526- 
66, 1900+14, 1627-41) were originally claimed to sit 
near but outside SNRs (ICline et ai~lll982t lAtteia et al.l 
[19871 IHurlev et al.lll999L I2000ir if the SNR and neu- 
tron star in these cases are physically associated and 
hence formed in the same supernova explosion, the large 
angular offsets of the SGR positions from the centres 
of the respective SNRs require high space velocities of 
thousands of km s . Indeed, it was hypothesized that 
magnetars received large kicks through magnetic field- 
drive n mechanisms that were ineffe ctive for ordinary pul- 
sars (|Duncan fe Thompson! Il992f) . However, in three 
of the four cases the density of SNRs on the sky in 
that region is sufficiently high that a chance projec- 
tion between SGR and SNR cannot be ruled out, while 
for SGR 1806—20, the classification of the extended 
nebula as a SNR is now th ought to h ave been erro- 
neous (|Gaensler et al. IfeOOlD . Recently, lEsposito et al.l 
( 2008) reported the detection of diffuse soft X-ray emis- 
sion around SGR 1627—41 and interpreted it as emis- 
sion from the supernova remnant, G337.00.1. Mean- 
while, direct measurements of magnetar proper mo 
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tions i mply much lower velocities than inferred b y these 
offsets (jHelfand et al. 1 120071: IRaplan et al.ll2009D . com- 
patible with the velocities observed for radio pulsars 
(I Arzoumanian. Chernoff fc Corde"s1 120021 iHobbs et al. I 
120051 : IChatterjee et al. Il2009f ). 

Assuming that the progenitor supernova explosion oc- 
curred at the geometric center of the SNR, the SGR 
has traveled approximately 80 arcminutes from its birth- 
place over its lifetime. For a distance to the SNR of 
800 ± 400 pc and an estimated age of 4000-7000 years 
([Leahy fc Tian 1120071 ) , the implied projected space veloc- 
ity for the SGR is 1300-6800 km s" 1 . As for the previous 
claimed SGR/SNR associations mentioned above, this is 
at or beyond the highest directly obs erved velocity ob- 
served for any other neutron star (cf. IChatterjee et al. I 
120051: iWinkler fc Petrel 120071 ). However, the proximity 
of this SNR allows a direct test: independent of dis- 
tance or projection angle, if the SNR and SGR are physi- 
cally associated, then the predicted proper motion of the 
SGR will be approximately 0.7-1.2 arcsec per year to 
the south, which should be measurable with the Chan- 
dra X-ray Observatory over even a relatively short time 
baseline. Confirming more associations would force a 
reconsideration of currently disfavored SGR-SNR asso- 
ciations, and of the highest velocities that neutron stars 
can attain. 
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